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ABSTRACT 

Three computer codes have been developed for the IBM-360 computer, 
namely: NUGAM2, NUGAM3 and SSALB. Two existing computer codes 
have been modified and made operational on the IBM-360 computer, 
namely: BREMRAD and OSS. The existing codes are described summarily 
and the new codes are detailed. 

The NUGAM2 code predicts forward and backward angular energy differen- 
tial and integrated distributions for gamma photons and fluorescent radiation 
emerging from finite laminar transport media. It determines buildup and 
albedo data for scientific research and engineering purposes; it also 
predicts the emission characteristics of finite radioisotope sources. It 
has been shown in this report that Sample results , in very good agreement 
with available published data, are presented. The code predicts data for 
many situations in which no published data is available in the energy 
range up to 5 MeV. 

The NUGAM3 code predicts the pulse height response of inorganic (Nal 
and Csl) scintillation detectors to gamma photons. Because it allows the 
scintillator to be clad and mounted on a 'photomultiplier' as in the experi- 
mental or industrial application, it is a more practical and thus useful 
code than others previously reported. Results are in excellent agreement 
with published Monte Carlo and experimental data in the energy range up 
to 4.5 MeV. 

The SSALB code determines single scatter albedo data analytically and is 
thus applicable to certain engineering problems. It was developed in 
support of the NUGAM2 code. 
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1 . . INTRODUCTION 

This final report was prepared by NUS CORPORATION for the National 
Aeronautics and Space Administration, Goddard Space Flight Center, under 
Contract NAS5-11781. The work scope consisted of modifying existing 
digital computer codes, namely; BREMRAD, NUALGAM and OSS. 

n 2) 

Under contract Task I, BREMRAD ' an existing IBM-7094 computer code 

for generating external and internal bremsstrahlung distributions, was 

modified. The modifications consisted of translation for operation on the 

NASA/GSFC IBM-360/91 computer, the redesign of code logic, the internal 

storage of frequently used data tabulations and the simplification of data' 

input. The completion of this task was reported to NASA/GSFC in Technical 

(3) 

Report NUS-TM-NA-109 and is thus not further discussed in this present 
report. 

(4) 

Under contract Task II, NUALGAM , an existing IBM-7094 computer code 

was modified. NUALGAM is a Monte Carlo code for calculating the emission 

of gamma photons from finite source media. The modifications resulted in 

the development of two new Monte Carlo gamma photon transport codes, 

NUGAM2 andNUGAM3, for operation on the NASA/GSFC IBM-360/91 

computer. This report is principally concerned with describing these new 

(5) 

codes. Also under Task II, the 05S computer code package was obtained 

from the Radiation Shielding Information Center, Oak Ridge National 
( 6 ) 

Laboratory' ' and made operational on the NASA/GSFC IBM-360/91. OSS is 
a Monte Carlo code for calculating fast neutron transport and pulse height 
response distributions in organic scintillators exposed to monoenergetic 
neutrons. This code is described in an appendix of this report. 
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Code NUALGAM was derived from an existing unpublished code — ALGAM, 

(7-9) 

developed at the Oak Ridge National Laboratories . ALGAM was 
developed for radiobiological application whereas NUALGAM was intended 
for the study of finite source photon number and energy distributions . 

Code NUALGAM predicts the angular-energy differential and integral escape 
of gamma photons from right-cylindrical homogeneous source media by 
means of the Monte Carlo technique. The modifications to NUALGAM 
reported herein consisted of the development of two gamma photon transport 
codes to predict buildup factors, albedos, fluorescent radiation emission 
and scintillation crystal pulse height response. The modified codes reported 
herein have been designated NUGAM2 and NUGAM3. 

The NUGAM2 code was developed to predict the forward and backscatter 
angular-energy distribution of photons and fluorescent radiation emerging 
from finite transport media of single or laminar composition. The code pre- 

9 • 

diets buildup and albedo data for either finite or semi-infinite transport 
media. An intermediate version of NUGAM2 , namely NUGAMl , was developed 
by NUS for NASA/GSFC as part of the S0SC code package^^^^. The shield 
transmission capabilities of NUGAMl have been extended and incorporated 
into NUGAM2 . 

The NUGAM3 code was developed to predict the pulse height response of 
inorganic scintillation crystals exposed to gamma photon radiation. It 
differs from previous codes which have been reported, in that the crystal 
model may be encased in a cladding material and mounted on a hypothetical 
photomultiplier . 

The present final report is composed of three volumes. Volume I describes 
the NUGAM2 and NUGAM3 Monte Carlo codes . Section 2 of Volume I 
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presents the theory upon which the code logistics are based while Section 3 
discusses significant sample results obtained. Section 4 consists of a 
summary. NUGAM2 and NUGAM3 are generally separated in Volume I 
through the use of subsection discussion, excepting where it is otherwise 
expedient. The important features of the OSS code as taken from reference 
(5) , are noted in Appendix I. 

Volume II of this report consists of a code users description for each of the 
codes NUGAM2 , NUGAM3 and OSS. The OSS code information is given in 
Appendix V of Volume II. In addition to a users description. Volume II 
also gives listings of code sample input and output; the code FORTRAN 
listings are omitted so that upon request a latest version may be supplied. 

Volume III of this report considers the statistics of the Monte Carlo method 
relative to the interpretation of the NUGAM2 and NUGAM3 results. A num- 
erical experiment using NUGAM2 code is presented and the results are 
interpreted. 



2 . 


THEORETICAL DISCUSSION 


2 . 1 General 

The NUGAM2 and NUGAM3 codes predict gamma photon transport distributions 
according to the Monte Carlo technique. Source photons are initiated stochas- 
tically and followed throughout their history in the transport medium. 

The NUGAM2 code categorizes those photons which escape from the medium. 

The escape categorization is performed as a function of escape energy and 
angle with additional separation into the backward and forward directions. 

Thus the code predicts "all space", albedo and buildup angular-energy 
distributions as a function of source photon energy and transport geometry. 

The distributions are predicted for photon number, energy and exposure. 

In addition to categorizing Compton and pair production photon escape 
the NUGAM2 code categorizes the escape from the transport medium 
of the characteristic x-rays of designated elements in the medium . 

The code also categorizes unscattered, single and multiple scatter escapes. 

The NUGAM3 code categorizes the energy-differential distribution of photon 
energy absorbed by a scintillation phosphor. It differs from previously 
published Monte Carlo scintillator response codes in that the phosphor 
may be encased in a cladding material and mounted on a "hypothetical" photo- 
multiplier tube. The code predicts the energy absorption spectrum and its 
Gaussian distributed counterpart analogous to the experimental measurement 
obtained with a multichannel pulse height analyzer. In addition to pulse height 
response the code predicts the angular-energy distribution of photons escaping 
from the "detector configuration" . The escape categorization is similar to 
that of the NUGAM2 code . 
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The NUGAM2 code Is designed for either a right-cylinder or a rectangular 
parallelepiped transport medium. The NUGAM3 geometry is right-cylindrical. 
The multiple media in NUGAM2 are parallel discs or slabs. The multiple 
media in NUGAM3 consist of a scintillation phosphor eg. Nal or Csl, clad 
on the front and sides and backed by a cylinder-simulated photomultiplier 
tube . 

The codes are designed to expose the transport medium to source photons 
originating from either a point isotropic source located outside or inside 
the medium (the internal source may be either a point or a uniform distri- 
bution) or parallel incident source photons. 

Since the codes do not at present follow photo, Compton or pair electrons 
and the bremsstrahlung they generate along their torturous transport paths, 
the codes are upper-energy limited. NUGAM2 is valid to ~ 6.0 MeV and 
NUGAM3 to ~4.5 MeV. Above these energies the codes will operate but 
will tend to deviate progressively from a "correct" prediction. 
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2.2 Monte Carlo Gamma Photon Transport Theory 

Codes NUGAM2 and NUGAM3 are programmed to predict gamma photon 
number transport in multiple media geometries. They employ the Monte Carlo 
technique of following and catagorizing a large number of photons from "birth 
to death"^^^' . They use random number and probability theory combined with 
known interaction distributions to determine such as source and collision 
site spatial location, as well as trajectory energy and direction throughout 
each history. For the purpose of clarity the text in this section assumes and 
thus refers to a right-cylindrical transport geometry, unless otherwise noted. 


After n Interaction events in the transport medium, photon number state 
may be characterized as 

= N (En/®n' ^n' ^n' ^n' ^m' ym' ^m^' 

where 


E 

e 

$ 

x,y,z 

n 

m 

o 


photon energy 

polar angle 

azimuth angle 

Cartesian co ordinates 

subscript to denote n^^ event 

n-1 , a subscript; 0, 1 . . . . , m, n, . . . . 

subscript to denote photon history origin, i.e. , 


which describes the spectrum at the arbitrary point Pj^ (x^/ yn< inter- 

action site Pn is outside the boundaries of the cylinder of length L and radius 
p, and Pjjj is within, then the "fate” of the photon is deemed as "escape", and is 
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so tallied by the code. If an escape is recorded at where = P©/ then 
the photon escape energy and direction is unaltered and identical with the 
initial or birth values ie. uncollided escape. A typical escape history is in- 
dicated in Figure 1 , where escape is shown between and P 3 . 

The codes employ Cartesian coordinates and direction cosines to define 
source and interaction positions and the intermediate trajectories . The 
codes are presently designed for external collimated point Isotropic sources 
of photons, parallel Incident photons and internal point or uniformly distributed 
sources. The determination of their various coordinates and direction cosines 
is described in Section 2.2. 

The path length between interactions, is a function of energy, material 
composition and geometry. The codes determine path length according to the 
Monte Carlo technique^^^^ as 

-tj + -62 + • ' •+ 1 C“logeRi“(Mi'f'l + M2'^2 ^ ^ 

where 

= the i^h random number in a sequence of random numbers, 

-tj = actual uncollided path length of photon in medium region j , 
between events m and n, 

Xj = total mean free path for a photon of energy E in medium j, 

Mj = total macroscopic interaction cross-section for a photon of energy 
E in medium j = Xj ^ . 

Figure 2 illustrates a typical photon trajectory between interactions in a 
a three-medium transport geometry. 
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For a single transport medium , j = 1, and thus Equation (2) reduces to 

-t = -XilogeRi . (3) 

Source photons are initiated with a probability of existence or weight W, 
equal to 1.0. This weight is reduced after each interaction by the ratio 
of the scattering to total cross-section. Thus photons are not lost to 
absorption unless their reducing weight drops below an assigned value, 
for example, 10“^. Photons maybe lost to absorption if their degraded 
energy drops below a "user" designated threshold value, eg. 0.025 MeV 
cutoff. 

Energy deposition within the transport medium at the n^^ interaction site is 
determined as 

Wm 

In =-^ C(Eni-l-02) CTpp+Em - En)] , (MeV) (4) 

where 

Wm 
Em 
En 

O’ 

pp,pe,s,t 
f^t 


= weight after m^^ interaction, 

= photon energy before interaction (MeV) , 

= photon energy after interaction (scattered photon), (MeV), 

= microscopic interaction cross-section, 

= subscripts to identify pair production, photoelectric , scattering 
and total cross-sections, 

= ^pp + ‘^pe + <^s • 
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Equation (4) neglects electron transport effects and is thus only totally 
valid where they are not significant eg. away from geometrical boundaries . 

The total macroscopic Compton scattering cross-section per atom is defined 
as 

“ ^i * *^es ' (5) 

where 


Zj^ = material or element atomic number , 

Ces = total microscopic Compton scattering cross-section per 
electron . 

The scattering cross-section cr is obtained from the Klein-Nishina re- 

6 S 

lationship as^^^^ 

= + , ( 6 ) 

where 

rQ = 2.8183 X lO'^^cm; classical electron radius, 

E = photon energy in moc^ (=0.51097 MeV) units. 

The macroscopic cross-section /it is determined from the microscopic 
cross-section cr in accord with the relationship 

^ = Ao , (7) 

i=l 
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where 


24 

Aq = Avogadro number; 0.6023 x 10 / 

p = material or element density , 

A = material or element atomic weight/ 

i = subscript to denote summation for elements ; i = 1 , 2 , . . . , k . 

Photon direction and energy after a scattering are governed by the angular- 
differential form of the Klein -Mishina distribution, with energy related to 
direction by the Compton scattering relationship^^ 

En = Em/d + (Em/0*51097)(l-coses)), (8) 

where 

©s r photon scattering angle, 

Em = photon energy before scatter (MeV) , 

En = photon energy after scatter (MeV) , 

The codes select scattering angle and energy at each interaction site in 
accord with the algorithm of Kahn outlined in reference (18) (reproduced 
in Appendix II ) . 

At photon energies greater than 2 moc (=1.02 MeV) the codes account for 
the pair production phenomenon. The codes temporarily store the parent photon 
characteristics and initiate a daughter photon with isotropically selected 
direction and energy mQC^ (=0.51097 MeV). Upon termination of the 
daughter history, the parent history is recontinued. In code NUGAM2 
the daughter history is initiated with twice the parent weight in order to 
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simulate an actual photon pair. In NUGAM3, the scintillation phosphor 
code, two back-to-back daughters are initiated, each with the parent 
weight, and followed separately to termination. 
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2.3 Geometry 


The NUGAM2 and NUGAM3 codes are designed for five different photon 
source geometries , namely: 

I. external axial point Isotropic source, 

II. external parallel Incidence, 

III. cases I and II with collimatlon, 

IV. Internal point Isotropic source, 

V . Internal uniformly distributed source . 

The various geometries allowed are Indicated In Figure 3 . The angle of 
parallel Incident photons may be user designated. 

For the case of an axial point Isotropic source as In Figure 3a, all 
source photons are Initiated at the coordinates (Xq, y^, z^) = (0,0, -r) 
with direction cosine 


= 1 - (1-r) Ri , (9) 

where 

T = r/(r^ + p2)l/2 ^ 

p = radius of transport medium cylinder , 

r = source-to-cyllnder distance, 
th ' 

= 1 random number in a sequence of random numbers . 

The azimuth angle is obtained as 

0o = (2Ri+i-l) . (10) 
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Direction cosines OLq and /3 q are determined as 

O^o = cos0o/q , (11) 

and 

pQ = (q^-r^)l/2 sin0Q/q , (12) 

where 

q = r/Vo . 

Collimation of the source cone is obtained by replacing the transport cylinder 
radius Pwlth p', the radius of collimation. 

Parallel incident photons are initiated on the front face of the transport 
cylinder, the z = Zq = 0 plane, at a radial distance d© from the axis and at 
a radial angle 0 q where 

1/2 

do = Ri • P> (13) 

and 

^o = (2Ri+i-l) . (14) 

If the parallel beam is collimated p replaces p. Equation (14) is identical 
to Equation (10). The directional cosines are (0£o, Pq, >b) = (0,0,l). 
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For rectangular-parallelepiped transport geometry exposed to an axial point 
source, the above equations are valid where pis the radius of the circle 
circumscribing the rectangular front face. In this geometry an additional 
test is required to ascertain if the individual photons "hit or miss" the 
rectangular front face. For the special case of axial incidence the 
coordinates (xo,yo#Zo) "user" specified. For uncollimated parallel 
incidence Zq = 0 and 


Xo 

w . Ri , 

(15) 

yo 

, ^ • Ri+ 1 / 

(16) 

where 



w = 

rectangular width, 


b 

rectangular breadth 


For the internally located point source case the photons are initiated at 
the specified coordinates (xQ,yQ,zo). For the internal uniformly distributed 
cylindrical source the coordinates are chosen by random number selection as 


Zq — R| . L , 

Xq = dQ COS0Q , 
Yo = ' 

(17) 

(18) 
(19) 

where 

■Jo = p • "li-f ' 

(20) 


00 ~ (2Ri+2"^) ^ ' 

(21) 
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and 


L = cylinder length . 

The initial direction cosines of each history for the internal sources are 
determined as 

2Ri -1 / 

cos[( 2 Ri+i -i)ir]. (l-ao^)^/2^ ^ 22 ) 

sin[( 2 Rj^j -l)ir]. (l-a^)l/2^ 

For the internal uniformly distributed rectangular source the coordinates 
{Xq,Vo,z^) are obtained from equations (15), (16) and (17). 

• 

When the photon path length I as determined from Equation (2), exceeds the 
vector path in the transport geometry the photon is deemed to have escaped 
from the medium geometry. In this event the history is terminated and 
categorized with respected to its energy, direction and weight. At the 
termination of each history in the NUGAM2 code, the energy and exposure 
(dose) is tallied. At this same point in the NUGAM3 code the absorbed 
energy is weight distributed and tallied. 

Although the code determines exact angles and location of photon escapes 
with respect to the transport medium boundary surface, a general categori- 
zation is carried out with respect to specific coordinate points, namely: 
the axial end points and geometric center. These locations (0,0,0), (0,0, L/2) 
and (0,0, L) are used as references points for albedo, 4ir and buildup data 
tabulations . 


^o = 
>o = 


Ik. 
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If a detection band about the Z-axis is assumed to be relatively distant 
from the transport medium then photons escaping from any point on the 
medium surface and striking the band, may be considered as having the 
same Z-axis directional cosine, ie. the medium is viewed as if a point source. 
If the detection band subtends the differential solid angle dO at (0, 0, O^z^l) 
then, the photon escape state may be characterized as 

N = N(E, fl), (23) 


where 


E = photon escape energy. 

In Equation (23) the understood arguments and subscripts are omitted. Multi- 
plying by differentials in energy and solid angle, the differential angular- 
energy photon number escape spectrum is given by 

dN (E f7.) = N (E, Q) dEdO , (24) 

keeping in mind that 

dfi = Sin0d 6d $, (25) 

where 


6 = polar angle between escape vector and Z-axis; 0^ 8^n , 

= azimuth angle of escape vector; 0 ^ 2 77. 
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Substitution of Equation (25) into Equation (24) and integration over energy 
and angle between desired limits, gives the total photon number escaping. 
The geometry for photon history escape classification is illustrated in 
Figure 4 for the reference point (0, 0, L/2). The escape classification 
geometry for the reference points (0, 0, 0) and (0, 0, L) is illustrated in 
Figure 5. Figure 5 gives the hemispherical segmentation for categorization 
of escape direction. Figure 6 gives the coded equivalent of Figure 5. 


17 



2.4 X-Ray Transport 

Code NUGAM2 is programmed to predict the escape from the transport medium 
of characteristic K x-rays originating in the photoelectric interactions. The 
present code version allows this option only for geometries consisting of a 
single transport medium. An additional restriction limits the x-rays to those 
of a single element in a medium composed of one or two elements . For example 
the escape of either the iodine or sodium K x-rays in a sodium-iodide medium 
may be predicted providing no adjacent media are considered. 


The number of element-B K-shell x-rays, generated at the n^h interaction site 
in a medium composed of elements A and B, is obtained as 


^P^B ^Kb 

^PPA+B ^P6B 


W, 


m 


B 


(2 6 ) 


where 


‘^Pee = 
^PSA+B " 

°kb 


microscopic photoelectric cross-section for photons of energy 
E and element B, 

microscopic photoelectric cross-section for photons of 
energy E and elements A+B, 

microscopic photoelectric cross-section for photons of 
energy E interacting with the K-shell electrons of element B, 
fluorescent yield of element B, 

weight of interacting photon prior to n^ interaction i.e. = 


The cross-section is obtained from reference (17) if the cross-section 
CTpgg is known. Table 1 reproduced from the reference, gives values of the 
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ratio O’pe/'^K ® function of element atomic number. The code determines 
element ratios by quadratic interpolation of the Table 1 data . Fluorescent 
yield data as a function of atomic number, reproduced from reference (19), 
are given in Figure 7 and Table 2. The code obtains Fg from the Table 2 
data by quadratic interpolation. 

The code initiates K series x-rays at the n^^ interaction site with a weight 

Wjjj, which is shared by the Ka2» in accord with relative 

intensities given in Table 3 (and input by the user). The code initiates K 

series x-rays in the same random isotropically selected direction with 

energies corresponding to the difference in the K and L or M shell energy 
( 20 - 22 ) 

levels . Intensity and energy data reproduced from reference (20) are 

given in Table 3. 

The code assigns a zero probability of x-ray escape for those interaction sites 
located at greater than a critical distance, d = 4.6X (Baj) from the transport 
medium bounding surface. For interaction sites located at less than the 
critical distance from the medium boundary, escape or absorption is 
determined by comparing the Monte Carlo path length obtained from Equation 
(2) with the vector distance to the boundary. 
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2.5 Scintillation Detector Response 

Code NUGAM3 determines the response of a scintillation crystal to incident 

\ 

gamma photons according to the Monte Carlo technique. The crystal may be 
encased in a cladding material and backed by a hypothetical "photomultiplier" 
analogous to the experimental situation. Thus the code allows the transport 
medium to consist of three distinct submedia, namely: 

I. front and side cladding 

II . " Photomultiplier" 

III. Scintillation crystal 

The code may be used for the combinations: I and III, II and III or III only. 

A typical detector geometry is shown in Figure 8 . 

The source photons are followed from their entry into the detector until they 
escape from the detector. The energy deposited at each photon interaction 
within the crystal site is weighted ,and subtallied as absorbed energy. The 
weighted subtallies are stored as weighted counts in the appropriate energy 
bins or channels at the termination of each history, i.e. after either escape 
from the detector system or weight/energy cut-off termination. 

A hypothetical photon history is shown schematically in Figure 9 . In this 
example the main path, indicated by the solid lines, is characterized by 
four interactions followed by escape from the detector. Photon weight is 
reduced at each interaction by the ratio of scatter to total cross-section; 
photon energy is also reduced at each interaction in accord with the Compton 
scattering relationships (see Section 2.2). Since the photon energy prior 
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to interactions 1 and 2 exceeds 2.0 mc^, the pair production threshold 
energy, two pair of "back to back". 1.0 mc^ photon histories are initiated 
as indicated by the dashed lines. In the figure the photon weights and 
energies (in parentheses) are shown before and after each collision. For 
example, it is noted that the main history initial energy and weight are 
8.0 mc^ and 1.0, respectively; the final (escape) values are 0.5 mc^ and 
0.1, respectively. 

The NUGAM3 code follows both main and annihilation histories; in the 
Figure 9 example five histories are followed. Figure 10 demonstrates 
the energy absorption and escape distribution corresponding to Figure 9. 

Figure 10(a) corresponds to the main history and Figures 10(b) and 10(c) 
correspond to the two pair of annihilation photon histories. Energy distribution 
matrices corresponding to Figure 10 are given in Table 4, 5 and 6. 

The energy absorption termination distribution for the example is given in 
Figure 11. The data for this figure, as tabulated in Table 7, is obtained 
from the "weights" and "absorptions per photon" given in Tables 4, 5 and 
6. The "absorption per photon" data of Table 5 is 6 mc^ larger in Table 7 
to account for pair electron energy deposition in the crystal. There is 
no electron energy addition for the second interaction since it occurs in 
the "photomultiplier" . Table 7 data of the same energy has been "grouped" 
for plotting in Figure 11. 

The NUGAM3 code tallies each history as in Figure 11. The result for a 
large number of histories is a spectrum of scintillation crystal absorbed 
energy. The solid points in Figure 12 present the Monte Carlo absorbed energy 
spectrum for a 3" x 3" Nal (T-t) crystal as determined by NUGAM3; the source 
energy in this example is 2 . 615 MeV. The solid points in Figure 12would be the final 
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useful result if experimental scintillators responded discretely to absorbed 
energy deposition. However, since the response is not unique but is 
instead spread out (smeared) over a range of pulse-heights, a "line -broadening" 
operation is necessary for simulation of the experimental spectrum. The line 
distribution in Figure 12 represents the "line -broadened" spectrum; it is 
discussed below. 

Experimentally it is an observed and fortunate fact that the "line-broaden- 
ing " in scintillators is approximately a normal , i.e. , Gaussian, distri- 
(23) 

bution . It has been found empirically that the standard deviation 
of the Gaussian, though not constant with energy, has a systematic depen- 
dence. In the NUGAM3 code the Gaussian standard deviation ct(E) is 
computed from the relationship^^^ ' . 

ct(E) = Ice"^ , (27) 


where 


k and n are determined by fitting the photopeaks of standard 
source spectra with Gaussians (k and n are the intercept and 
slope in a log-lot plot of c;(E)) . 


The NUGAM3 code redistributes the j = 1 to m channel elements of the 
Monte Carlo absorbed energy spectrum M(Ej) , over i = 1 to n channels 
according to the Gaussian relationship (defined in reference (26)) 


9ij 


(27T)~^ M(Ej) 
<7(Ej) 


exp 




1/2 


Ej - Ei 


ty(Ej) 


) • 


( 28 ) 


The final Gaussian broadened spectrum is obtained from the summation 
of the m spectra as 
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(29) 


m 

Gi = ^ gij ; i = l,n. 

J=1 

‘In equations (28) and (29) and are the energies at the channel mid- 
points . In the event that the standard deviation o(E) is small relative 
to the selected constant channel width, Av = Em/m, the code determines 
the Gaussian in each i^^ channel as a histogram i.e. Equation (28) is 
integrated over each channel from V-0. 5 Av to V+0.5 Av^^^\ The 
distribution in Figure 12 is the Gaussian smoothed absorbed energy spectrum 
corresponding to the solid points. 


In addition to determining the Monte Carlo absorption and Gaussian smeared 
spectra the code also calculates the detector photofraction^^®^, P(E) 
and efficiency, e(E), according to the relationships 


P(E) 


Full Energy Absorption Counts 
Total Counts 


(30) 


and 


_ Total Counts . (31) 

Number of Detector Incident Photons 

The "number of detector incident photons" is defined by the frontal cladding, 
i.e. , a lesser number are defined by the crystal front face for obvious 
geometrical reasons. 
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3. DISCUSSION OF RESULTS 


3 . 1 Introduction 

In this report section, results obtained with the NUGAM2 and NUGAM3 codes 
are presented and discussed. The NUGAM2 code was used to obtain gamma 
photon buildup, albedo, finite source self-absorption and finite source emission 
data. The NUGAM3 code was used to determine Nal(Tt) scintillation detector res- 
ponse function data. The NUGAM3 code was also used to obtain encapsulated 
source self-absorption and emission data. 

Buildup factor and albedo distributions were determined for gamma photons 
incident on finite laminar shields using the Monte Carlo technique. Angular 
energy differential and integrated number, energy and exposure distributions 
were obtained as a function of angle and energy (0.5 to 6.0 MeV) of incidence 
as well as dimensions of the medium. The buildup factors for single medium 
composition are in good agreement with the published theoretical^ and 

(31) 

experimental data . The albedo data for single media are in similar good 
agreement with published data^^^'^^K The albedo data is presented in 
Section 3.2. The buildup data is presented in Section 3.3. The source 
self-absorption and emission data is given in Section 3.4. 

NaUT-f-) experimental pulse-height response functions were numerically 
simulated using the NUGAM3 code. The analysis was carried out for 
right-cylindrical Nal crystals encased in typical side and front cladding 
material and rear-mounted on a material typifying a photomultiplier assembly. 

The results are for exposure of the encased detector to isotropic point sources 
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of gamma photons in the energy range up to 4.45 MeV; parallel beam photon 
incidence was also investigated. The results obtained were: crystal energy 
deposition and Gaussian-smeared pulse-height response distributions; 
photo fractions; interaction efficiencies and detector angular-energy differ- 
ential escape distributions. The importance of single and double scatter 
effects was also investigated. The results for a "bare" crystal are in ex- 
cellent agreement with the data of reference (44). The results for the clad 
crystal simulate the experimental data of reference (23) with very good agree- 
ment. 


3.2 Albedo Results 

The NUGAM2 code determines gamma photon albedo data as a function of photon 
angle of incidence and energy, material atomic number and composition, 
laminar disposition, material lateral dimension and thickness, and incident 
radiation beam geometry such as parallel beam , point source , etc . 

The number of gamma photons backscattered as a function of the transport 

medium total thickness L and within an escape category radius d located on 

the transport medium backscatter plane, as in Figure 13, were determined to 

ascertain the limitations of the albedo technique . Total number albedos for 

0.662 and 1.25 MeV normal photon incidence are given in Figures 14 and 15 

as a function of both L (open points) and d (solid points) , in units of mean 

free path X (E) , with transport medium atomic number Z as a parameter. The 

arrows on the right hand ordinate of the figures indicate comparison data from 

the references (33-43). For Z = 13 the number albedo reaches a saturation 

value at d = d„ < 0. 25 X for both 0.662 and 1.25 MeV. As the figures indicate, 

approximately half of the backscattered photons emerge within a characteristic 

distance, d = dj^, typical values of which are given in Table 8 in units of 
2 
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X (E). Analogous half and saturation thicknesses Lj^ and Leo were observed 

2 

to be approximately 85% of the dj_ and doo values. The total number albedos 

2 

may be approximated as a function of either d or L, by the exponential 
relationship form reported by others. (45, 46) The saturation total number, 
energy and exposure albedos were found to be in excellent agreement with 
the published Monte Carlo data and in good agreement with the experimental 
data. 

Figure 16 gives total number and exposure albedo distributions for normally 
incident 1.0 MeV photons as a function of lamination thickness and composi- 
tion, for L = iX, 2X, and 4X. The number (open points) and exposure (solid 
points) albedo are seen to vary over the entire lamination thickness fraction 
range for A1 - Pb and A1 - Fe . The lamination order in the figure is source - 
A1 - Fe (or Pb) . For example , the total number albedo for 1.0 MeV photons 
is 0.08 for iX aluminum 'backed' by 3X of Pb (thickness fraction = 0.25/0.75); 
iX A1 = 6.034 cm, 3X Pb = 3.865 cm, L = 4X = 9.899 cm. 

Figure 17 gives the total number albedo distributions for normally 'incident' 

1.0 MeV photons and a medium thickness L = 4X as a function of lamination 
arrangeifient order, e.g. , A1 - Pb and Pb - Al. The albedo example in the 
last paragraph reduces to 0.01 if iX Pb is 'backed' by 3X Al. In fact it can 
be seen that the albedo is that of Pb only until the Pb thickness is reduced 
to less than 0.25X. The same behavior is evident for the Al - Fe and Fe - Al 
data. 

Figure 18 shows the code calculated total number albedo as a function of 
normally incident photon energy, with atomic number as the parameter. 
Calculations were not carried out above 6.0 MeV since the code does not 
account for bremsstrahlung generation which will be significant at higher 
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energies. Similarly the Pb data In the figure does not account for fluorescent 
radiation backward emission; this Is significant for high Z as indicated in the 
references (33, 34). 

Figure 19 gives examples of the total number albedo as a function of the angle 
of incidence 0 o; do is measured with respect to the outward normal of the 
scattering medium. For the sake of the expected albedo at 6 o = 90°, the 
albedos in this figure are multiplied by cos 6 o. The increasing albedo with 
increasing angle 0o is very evident in the figure for Z = 82, even though it 
does occur for Z = 13. The dependence on energy is also noticeable. 

Figure 20 shows the ratio of the 'single scatter' to 'total scatter' for the 

total number albedo, as a function of energy. The ratio oc is 

Nss NT 

plotted as a function of E^, with atomic number as a parameter. Typical 

single and double scatter to total number albedo ratios oc and a , are 

ss ds 

given in Table 9 . It can be seen that multiple (>2) scattering in lead accounts 
for only £ 1% of the number albedo at E ^ 0.662 MeV. However, even for 0.412 
MeV photons incident on aluminum, multiple scatter accounts for ~ 50% or 
more of the total number albedo. The Monte Carlo single scattering data 
were substantiated by a double numerical integration of the analytic differential 
transport equation, indicating that albedos for high Z media may be readily 
approximated by this method. The analytic transport equation used in this 
evaluation was defined, with reference to Figure 21, as 

^ax Q 

a( 0 ,Eo)e-^s^^^°s sin 6 de dz, (32) 

8=0 
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where 


r 


e 


Ne 

o-(6,Eo) 

Ms 

Eq 

Es 


classical electron radius, 

2.8183 X 10"3 (cm),^^®^ 

3 

number of electrons/cm , 

(2 6) 

Kleln-Nishina Unction (cm^/ster. -electron), 

total cross-section for E© (cm“^), 

total cross-section for Eg (cm"^), 

incident photon energy (MeV), 

; ©s = - ®/ (MeV) . 

1 + (Eo/0.511)(1-cos6s) 


Table 10 presents a comparison of Al, Fe, Sn and Pb analytic (Equation (32)) 
and NUGAM2 Monte Carlo single scatter albedos for different values of 
thickness L, escape category radius d and energy E (0.662 and 1.25 MeV). 
The agreement can be seen to be very good for low atomic number and as poor 
as 15% "error" for high atomic number. The deviation expresses to a large 
extent, the statistical error in the Monte Carlo data. The analytic equation 
was encoded and solved using Simpson's Rule in a separate computer code: 
SSALB. Code SSALB is detailed and listed in Appendix VI of Volume II of 
this report. 


Code SSALB is considered as the basis of a very useful engineering code for 

radiation backscatter analysis . For example. Table 10 indicates (ajyj ) 

s s 

anal. = 0.0656 for Z = 13 and E = 0.662. Table 9 gives R = 0.24, thus 

o ss 

using the analytic single scatter albedo the total albedo may be estimated 
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as = (0.0656/0.24) = 0.274, which is in very good agreement with 0.26 
in Figure 14. This approach based on Monte Carlo data, offers an economical 
method for engineering analysis . 

3.3 Buildup Results 

The NUGAM2 code determines gamma photon buildup distributions as a 
function of photon angle of incidence and energy, material atomic number 
and composition, laminar disposition, material lateral dimension and thickness, 
and incident radiation beam geometry. 

The energy penetration buildup of 1.0 MeV photons incident on aluminum, 
iron and lead media, as a function thickness, is compared directly in Figure 22 
with the data of reference (29) and indirectly with reference (30). The 
reference (29) data which in turn compares well with the reference (30) data 
is indicated by the solid line curves while the NUGAM2 data is given by the 
points. The reference (29) data was determined by the Monte Carlo method 
as applied to semi-infinite media. The NUGAM2 data was determined for 
finite media of "large" lateral extent. The agreement is seen to be very good. 
NUGAM2 data for aluminum is compared with reference data for ordinary 
concrete . 

For normal axial incidence it was observed that ~ 99% of the transmitted 
scattered photons emerged from the transport medium within a characteristic 
radial distance doo approximately equal to the shield thickness. Figure 23 
gives photon number penetration buildup factors for L = 0.5X, iX, 2X, 4X 
and 6X as function of radial distance d (in units of X(E)) from the axis of a 
cylindrical tin shield axially exposed to 0.5 (open points) and 2.0 (solid 
points) MeV photons. The figure indicates that photon emergence from a 
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shield is 'negligible' at perpendicular distances >3X from a normally incident 
axial photon beam. Figure 24 gives photon number penetration buildup 
factors for L = 2X and 4X as a function of d for H 2 O, Al, Fe, Sn and Pb 
exposed to 2.0 MeV photons. The behavior is seen to be similar to that 
given in Figure 23 for tin. 

Number and energy buildup factors for photons of 1.0 MeV, normally incident, 
on shields of thickness L = 4X and 6X, are given in Figure 25 as a function of 
lamination thickness for Al - Pb, Pb - Al and Al ^ Fe. The distributions in 
this figure are in agreement with qualitative arguments based on the photo- 
electric effect and Compton scattering. The energy buildup factors are less 
sensitive than the number buildup factors to the shield geometry and composi- 
tion, especially when the media are relatively proximate in atomic number 
as in the case of Al - Fe. The Al - Pb data indicates that a laminar shield 
of 3X aluminum, backed by iX lead, has energy and number buildup approxi- 
mately equivalent to a shield of 100% lead. 

Figure 26 gives the energy spectrum of scattered photons penetrating Pb - Al 
shields of thickness L = 4X as a function of the lamination fraction for 2.0 MeV 
incident photons; the unscattered penetrating photons are not shown. The 
all Al shield exhibits a low energy 'buildup' whereas the all Pb shield exhibits 
a low energy 'depletion,' which is in agreement with qualitative considerations. 
The effect of reverse lamination order is seen in Figure 27 where the spectra 
of Pb - Al and Al - Pb are shown; shield total thickness and lamination 
thickness fractions are the same. Figure 27 shows that the buildup is a 
minimum when low Z is 'backed* by high Z material; the high Z acts as a 
'catcher' for low Z scattered photons. The spectra of Figures 26 and 27 are 
consistent with experimental spectra of reference (31). 
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3.4 Source Self-Absorption and Emission 

The NUGAM2 code determines gamma photon self-absorption and emission 
from either cylindrical or rectangular finite source geometries . Initial 
isotropic source distributions may be specified either as either 'at a point' 
or uniformly distributed. The present version of NUGAM2 allows the source 
to be in one lamination of the laminar geometry, and thus unclad "sidewise" . 
A NUGAM3 option allows the substitution of a source medium for the NaKT-t) 
crystal medium and further allows the substitution of source encasement for 
the normally used crystal cladding and photomultiplier. Thus, NUGAM3 
may be used for encased source predictions. 


Figure 2 Spresents predicted unscattered and scattered photon number fractions 
emitted into 4n steradians for a 2 x 2 in. samarium-oxide cylindrical source 
of density 1.5 gm/cm^, as a function of source energy. The figure also gives 
the unscattered fraction for a 1 x 1 in. source; the scattered fraction is omitted 
since it lies close to that given for the 2x2 in. source. The effect of the 
samarium K-electron absorption edge at 46.8 keV is evident for both sources. 

Figure 2 9 presents unscattered and scattered photon number fractions for 
escape into 4ff steradians as a function of source size for selected source 
energies and materials. For a 0.5 MeV photon source homogeneously distri- 
buted within an 8 X 8 cm aluminum cylinder, it is seen that the fraction of 
scattered escapes approaches that of unscattered escapes. The figure also 
shows that the scatterered fraction for aluminum is approaching a maximum. 

This maximum, followed by a gradual decrease, was observed at ~ 4 x 4 cm 
for the samarium example and is related to the 'quenching' effect of photo- 
electric absorption, which is increasingly pronounced with increasing atomic 
number and density. During the course of this work, the differential angular 
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and energy variation of photon escape was obtained as a function of atomic 
number, transport medium density and cylinder helght-to diameter ratio. 
Integral data was determined from the differential results. 

3.5 Scintillation Crystal Results 

The NUGAM3 code determines the pulse-height response distributions of 
a scintillation crystal exposed to a point source or a parallel beam of 
gamma photons. The code optionally allows the crystal to be clad and/or 
'photomultiplier' mounted. 

Table 11 gives photofractions and detection efficiencies determined by the 
NUGAM3 code for an unclad (u) (bare) and clad (c) (front and side cladding, 
and photomultiplier) 3" x 3" right cylindrical Nal(T't') scintillation detector 
exposed to parallel incident photons in the energy range 0.32 to 2,754MeV; 
a 4" X 4" crystal is included for data comparison at 4.45 MeV. 

The calculations assumed 0.125" front and side cladding and a hypothetical 

"photomultiplier" composed of aluminum with a density of 2.0 g/cc; a 

1.18 g/cm beryllium absorber was also assumed. The cladding dimensions 

(47) 

and composition were taken as those of typical commercial detectors 
Unclad Monte Carlo photofraction and efficiency data from references (44) 
and (48) are given in parentheses; clad Monte Carlo data has not been 
reported. The NUGAM3 and reference data can be seen to be in excellent 
agreement. 

Figure 30 compares the NUGAM3 code Monte Carlo generated spectral re- 
sponse of a 3" X 3" Nal(T'I') scintillation detector response to 10 cm distant 
51 85 22 

point sources of Cr, Sr and Na, with the actual experimental response 
reported in reference (23). The solid-points are the Monte Carlo data and 
the line distributions are as obtained by graphically connecting the digital 
data of the reference. The Monte Carlo and experimental spectra were 
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photopeak area and pulse-height normalized; pulse-height normalization 

was carried out by using the algorithm of reference (49). The reference (49) 

algorithm was also used to smoothen the Monte Carlo absorbed energy spectra 

prior to Gaussian smearing. The Monte Carlo data was 'smeared' using 

k arid n values in Equation (27) based on a fit of reference (23) photopeaks. 

(47) 

The NUGAM3 determination assumed 0.125" front and side cladding 

and a hypothetical "photomultiplier" composed of aluminum with a density 

of 2.0 g/cc; the front cladding was faced with a 1.18 g/cm beryllium absorber. 

As Figure 30 shows , the agreement between calculation and experiment is 
85 22 

excellent for Sr; very good for Na except in the 1.275 MeV photo peak valley; 
51 22 

and good for Cr. The Na valley discrepancy, although not large, may 
be removed by adjustment of the NUGAM3 'smearing' standard deviation ff(E) , 
given in Equation (27). The deviations in the ^^Cr Compton continuum are 
considered to be small and partly due to the statistics of the Monte Carlo 
determination. It is noted that other than normalization and smoothening 
no unusual effort was made to fit the calculations to the experimental data. 
Additionally some details of the experiment are not positively known. 

Figure 31 shows the Gaussian 'smeared' effects of front and side cladding, 
and 'photomultiplier', as determined by the NUGAM3 code for E = 0.6616 MeV. 
Response A consists of the backscatter distribution for the photomultiplier. 
Responses B and C consist of side and front cladding distributions. The sum 
distribution of A + B + C is also given. Responses for A, B, and C correspond 
to the geometry of Figure 30. In order to determine if the trends of the dis- 
tributions were real and since they were obtained from the difference of two 
distributions, i.e. , 'with and without A or B or C, ' an abnormally large 
thickness of cladding was investigated, namely: D and E in Figure 31. 
Responses D and E tend to confirm the general shapes of B and C. Additional 
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photomultiplier dimension introduced only a negligible change in Response A, 
indicating that an 'infinite' 'photomultiplier' had been assumed. 

Although the distributions of Figure 31 are for 0.6616 MeV photons the same 
general behavior was observed for the energy range 0.279 to 2.754 MeV. 

It was deduced that the photomultiplier backscatter effect is primarily an 
albedo spectrum which exhibits a 'single' scatter peak at an energy 
corresponding approximately to a photon scattering angle of ~ 150° to 170°. 
The side cladding effect is seen to be similarly an 'albedo' spectrum but 
with a broader and higher energy 'single' scattering peak since allowable 
scatter angles range from 0° to 180°. The frontal cladding effect is seen 
to be a distribution extending a/er almost the entire spectral energy range. 

Its contribution to the photo peak valley is apparent. 

Figure 32 presents the clad and unclad (dashed distribution) spectral 
response of a clad and unclad 3" x 3" Na^T-t) detector exposed to a 10 cm 
distant point source of 2.615 MeV gamma photons. The single and double 
escape peaks are seen in both spectra but the 0.51 MeV annihillatlon 
distribution is only evident in the 'clad' spectrum. This is one very important 
feature of the NUGAM3 code over previous codes, namely that it calculates 
the 0.51 MeV distribution. In 'thickly' clad scintillations crystals such as 
are proposed for space flights this is of considerable importance. 

Figure 33 presents the Monte Carlo calculated high energy portion of the 
Compton 'tail' on which the Gaussian photopeak 'sits,' as a function of 
source photon energy. The distributions are normalized to 50,000 source 
photons emitted in the solid angle subtended by the cladding of a 3" x 3" 
crystal; they are also pulse-height normalized to channel 100. Typical 
photopeak widths (at half-height) for 0.32 and 2.615 MeV would be ~ 10.0 
and ~ 4.6 channels on the normalized scale of Figure 33, respectively. 
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Figure 34 presents photofractions as a function of side and front cladding 
thickness as determined by NUGAM3 for a 3" x 3" NalCT-t) crystal exposed 
to a 10 cm distant point source of 0.662 MeV gamma photons. The behavior 
is consistent with spectra of Figure 31. The front cladding is seen to have 
a greater influence on photofraction than side cladding. 

Figure 35 presents relative photofraction determined by NUGAM3 as a function 
of the collimation of a parallel beam incident on a clad 3" x 3" scintillation 
detector, for E = 0.662 and 1.25 MeV. Collimation is expressed in a 
ratio: (radius of collimation)' (radius of crystal) . Thus the ratio range 
is greater than 1.0 if the cladding is taken into account. The energy 
dependence is evident from the separation of the two distributions in the 
figure. The absolute photofractions corresponding to the relative photo- 
fractions at the axis and cladding surface are given. The distributions are in 
very good agreement with the data given in reference (50). However, it is 
noted that actual collimation devices will influence photcfractions to the 
extent that photons scattered in the collimator intercept the detector. 
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4. SUMMARY AND CONCLUSIONS 


Five FORTRAN language digital computer code packages have been made 
operational on the NASA/GSFC IBM-360/91 digital computer. The codes 
are BREMRAD, NUGAM2, NUGAM3, SSALB, and OSS. The existing BREMRAD 
code was modified as described in the reference (3). The OSS code was 
obtained from the RSIC/ ORNL and made operational at GSFC; it is described 
in the appendices to this report and the references (5, 6). The NUGAM2 and 
NUGAM3 codes were developed from the NUALGAM code, to which they now 
bear very little resemblance. The SSALB code was designed and developed 
to substantiate the NUGAM2 code. Because of their complexity and 
originality the largest part of the work scope effort was necessarily devoted 
to the NUGAM2 and NUGAM3 codes. The remainder of this section refers to 
the NUGAM2 and NUGAM3 codes. 

The NUGAM2 code was developed to predict forward and backward angular 
energy differential and integrated distributions for gamma photons and 
fluorescent radiation' emerging from finite laminar transport media. The 
code is thus designed to determine buildup and albedo data for scientific 
research and engineering purposes. The code is also designed to predict 
the emission characteristics of finite radioisotope sources. Asymptotically 
the code will predict data for 'large' (~ semi-infinite) transport media. It 
has been shown in this report that sample predicted data is in very good 
agreement with available published data. However, since the code predicts 
data for many situations in which no published data is available , it is 
concluded that it is a contribution to the field of radiation transport. 

The NUGAM3 code was developed to predict the pulse height response of 
Inorganic scintillation detectors to incident gamma photons . Because the 


36 



code is designed to allow the scintillator to be clad and mounted on a 
'photomultiplier' as in the experimental or industrial application, it is a 
more practical and thus useful code than others previously reported. It 
is shown in this report that the code's results are in excellent agreement with 
available published data. However, since the code predicts data for 
application in those situations where reported data is not available , it 
is concluded that it is a contribution to the field of radiation spectrometry. 

Insofar as the NUGAM2 code is upper energy limited to the region of 5 to 
6 MeV, it is proposed that its capability be extended to include bremsstrahlung 
generation resulting from electron production in the transport medium. 

Since the code does not employ deep penetration biasing techniques it is 
proposed that this be considered to allow its application to thick (greater 
than 6 to 8 mean free paths) transport media. Finally it is proposed that 
the code be used to develop engineering and scientific basic data for 
terrestrial and space application. 

Although the NUGAM3 code was developed to encompass the radioisotope 
general energy range, up to ~ 4.5 MeV, it is proposed that it be extended 
to be useful up to ~ 10 MeV. This would allow it to be of considerable 
value in gamma photon astronomy and planetary spectrometry as well as for 
many terrestrial applications. It could be applied to determine complex 
detector responses where standard sources cannot be used. The proposed 
energy extensions require that the code be designed to consider electron 
escape and bremsstrahlung generation. Finally, it is proposed that the 
NUGAM3 code be used to develop a basic data library of response function 
distributions for scientific and industrial applications. 
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FIGURE 2 

TYPICAL PHOTON TRAJECTORY BETWEEN INTERACTIONS 
IN A LAMINAR GEOMETRY 
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FIGURE 4 

GEOMETRY OF PHOTON HISTORY ESCAPE CLASSIFICATION 
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FIGURE 5 

SEGMENTATION OF HEMISPHERE FOR 

DEFINING DIRECTION OF FORWARD AND BACKWARD SCATTERING 
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FIGURE 7 

GRAPHICAL SUMMARY OF K-FLUORESCENCE YIELDS*'*^* 
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FRONT CLADDING 


"PHOTOMULTIPLIER 



FIGURE 8 

TYPICAL SCINTILLATION DETECTOR ARRANGEMENT 
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ESCAPE PHOTON 

FIGURE 9 

HYPOTHETICAL PHOTON HISTORY FOR 
CLAD AND "photomultiplier"- BACKED SCINTILLATION CRYSTAL 



FIGURE 10 

ENERGY DISTRIBUTION OF FIGURE 9 
HYPOTHETICAL PHOTON HISTORY 
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INTERACTION NUMBER 



FIGURE 11 

PHOTON ENERGY ABSORPTION TERMINATION DISTRIBUTION 
CORRESPONDING TO FIGURE 9 HISTORY AND TABLE 7 DATA 



COMPARISON OF MONTE CARLO ABSORBED ENERGY SPECTRUM 
AND GAUSSIAN "SMEARED" SPECTRUM FOR 2.615 MeV PHOTONS AND 
3" X 3" CLAD Nal(TI) CRYSTAL 



CHANNELS 





FIGURE 13 
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TOTAL NUMBER ALBEDO AS A FUNCTION OF L AND d 
FOR Eq = 0.662 MeV AND NORMAL INCIDENCE 
FIGURE 14 
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TOTAL NUMBER ALBEDO 



FIGURE 16 

NUMBER(dashed) and EXPOSURE(solid) ALBEDOS for 1.0 MeV PHOTONS 

AS A FUNCTION OF LAMINATION THICKNESS FRACTION:. 
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TOTAL EXPOSURE ALBEDO 




TOTAL NUMBER ALBEDO 



FIGURE 17 


NUMBER ALBEDOS FOR 1.0 MeV PHOTONS 
AS A FUNCTION OF LAMINATION 
THICKNESS FRACTION 
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INCIDENT ANGLE, Go (Degrees) 

(TOTAL NUMBER ALBEDOO(|m)* COS ©6 AS A FUNCTION 
OF INCIDENT ANGLE e„. FOR Eo=0.662 AND 1.25 MeV 
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TIO AS A FUNCTION OF INCIDENT ENERGY E, 




GEOMETRY FOR ANALYTIC , 32 , 
SINGLE BACKSCATTER EQUATION 

FIGURE 21 
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ENERGY BUILDUP FACTOR 



SHIELD THICKNESS 
(MEAN FREE PATHS, X(E)) 

FIGURE 22 

COMPARISON OF PRESENT WORK 1.0 MeV ENERGY 
BUILDUP FACTORS FOR Al { + ), Fe ( • ) 

AND Pb { X ) WITH THE DATA OF PIEPER AND BEACHrJ^ND PEEBLEsl”’ 
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ESCAPE RADIUS, d (UNITS OF ME)) 

FIGURE 23 

NUMBER BUILDUP FACTORS AS A FUNCTION OF RADIAL 
ESCAPE DISTANCE d AND SHIELD THICKNESS L FOR 
0.5 (dotted) AND 1.0 MeV (solid) PHOTONS INCIDENT ON TIN 




NUMBER BUILDUP FACTORS 


6.0 



FIGURE 24 

NUMBER BUILDUP FACTORS AS A FUNCTION 
OF RADIAL ESCAPE DISTANCE d FOR 2.0 MeV PHOTONS 
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FIGURE 25 

NUMBER(dashed) and ENERGY(solid) BUILDUP FACTORS for 1.0 MeV PHOTONS 
AS A FUNCTION OF LAMINATION THICKNESS FRACTION 
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NUMBER OF TRANSMITTED SCATTERED PHOTONS 
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NOTE: (A) UNSCATTERED PHOTONS 
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FIGURE 26 


SCATTER SPECTRA FOR (2.0 MeV) PHOTONS TRANSMITTED 
BY Pb - Al SHIELDS AS A FUNCTION OF LAMINATION 
THICKNESS FRACTION 
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NUMBER OF TRANSMITTED SCATTERED PHOTONS 



SCATTERED PHOTON ENERGY, E(MeV) 


FIGURE 27 


SCATTER SPECTRA FOR <2.0 MeV) PHOTONS TRANSMITTED 
BY A1 - Pb AND Pb - Al SHIELDS 




NOI. 



0.01 0.1 1.0 10.0 

PHOTON SOURCE ENERGY (MeV) 

RGURE 28 

PHOTON NUMBER FRACTION EMITTED INTO 
STERADIANS FROM 51x1203 SOURCE CYLINDER 
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ESCAPE FRACTION 



SOURCE DIAMETER AND HEIGHT (CM) 


FIGURE 29 


PHOTON NUMBER FRACTION EMITTED INTO An 
StERADIANS AS A FUNCTION OF SOURCE SIZE FOR 
SELECTED INITIAL PHOTON ENERGIES AND SOURCE MEDIA 
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CHANNELS 


RGURE 30 

COMPARISON OF REFERENCE (23) EXPERIMENTAL SPECTRA 
AND WITH NUGAM3 MONTE CARLO SPECTRA 
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RELATIVE COUNTS / CHANNEL 





UNIT AREA 2.615 MeV 



70 


DISTRIBUTION FOR CLAD AND UNCLAD 3" x 3" Nal (T;Z. ) 
CRYSTAL EXPOSED TO POINT SOURCE AT lOctn DISTANCE 



100 



CHANNELS 
FIGURE 33 

COMPTON TAIL DISTRIBUTIONS FOR 3" x 3" 

Nal (7Z) CLAD CRYSTAL AND 10cm DISTANT POINT SOURCE 
(GAUSSIAN SMEARED AND NORMALIZED TO 50,000 PHOTON HISTORIES) 
(1.18 g/cm2 ABSORBER) 
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RELATIVE PHOTOFRACTION, P T 



COLLIMATION RATIO, b/R 


RGURE 36 

RELATIVE PHOTOFRACTION AS A FUNCTION OF 
COLLIMATION OF BROAD PARALLEL BEAM FOR 
CLAD 3" X 3'73‘gal (TjL) CRYSTAL 
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TABLE 1 


THE EXPERIMENTAL PATIO (T /tJ^ AS DERIVED 

pe K K-edge 

FROM KIRCHNER [ 1930 3 BY DAVISSON [ l965]/^^^ 


Element 

K-edge 

Element 

(T /t ) 

^ pe^ rK-edge 

i« 

1.00 

20^^ 

1.11 

4®" 

1.03 


1.12 


1.04 

29^'^ 

1.13 

7“ 

1.05 

42 Mo 

1.16 

8° 

1.06 

50®" 

1.18 

jlNa 

1.075 

53^ 

1.19 

12 '^® 

1.08 

74 '^ 

1.20 

13^^ 

1.09 

79"^ 

1.22 

14®‘ 

1.09 

32'’‘= 

1.22 

18^ 

1.10 

92 '^ 

1.235 

19^ 

1.10 
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TABLE 2 


K FLUORESCENT YIELD AS A FUNCTION OF ATOMIC NUMBER^^^^ 


Atomic Number 

Yield 

10 

.017 

15 

.050 

20 

.160 

25 

.30 

30 

.475 

40 

.71 

50 

.83 

60 

.90 

70 

.95 

80 

.967 

90 

.975 
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TABLE 3 

X-RAY CRITICAL ABSORPTION AND EMISSION 
ENERGIES IN KEV < 20 ) 


K S»rlM 


Um 

a 


®2 

»i 


'2 

\ 

S 

- 

Tr«a«UloB 
K ■. ■■ 


Sn 

Si 

Si. 

■Si 


■lA 

VS 


A^roila*i« Ist«M>Vjr 
(MUtlv* K^) 

150 

100 

50 

15 

5 

1 

i 

i 


L 


Um 

a 

“1 

“2 


»2 

») 






7 * 

'6 

PraMlttoB L. 


s 

V 

V 

s 

“m 

"ri 

®rv,if 


"r. 

*777 

“m 

V 

*ppraita»t« iBtBBBlt/ 
(RBUtlM I^) 

110 

100 

10 

80 

60 

p 

20 

P 

10 

10 

10 

10 

20 



■s. 



S 

s 

4 »b 

1 ^ 1 

*iXBb i \ Sg \ \ 


1-1 

i-u 

ti 

5-1 

o.oip 

0.0216 

0.055 

0 .U 6 

0.192 



0 . 

0 . 

0 . 

052 

uo 

105 


i 


6 ^ 

7 - 8 

8 - 0 
9 -r 

tO-M 

U-R* 

0.263 

0.599 

0.551 

0.607 


0.398 

0.523 

o. 6 n 

Eo^r— °:S'i -1 

i ! ‘ 

^Oi.O*8 0.022 ► 0.02? 

0 . 055 ~’ 0.05* ' oroj* 


IVAl 

iWsi 

U 5 -P 

IT-Cl 

Ifi-A 

19-1 

ao-c«_ 

fil^ae 

e2-Ti 


S 6 *M 

7 l<o 

a^tt 

29 ”Cu 

I 

ja-c* 

5 >H 

3 fc* 8 « 

5 *«>L 

36-lr 


1.30J ! 
1.559 ‘ 
, l.«38 ! 

*^ 2.^70 




i:S 

«- 5 Jr_ . . 

7.111 

7.709 

fi.jji ' B.jse 
6.960 e.976 
9 '^ 

10. ^ ■ 10.J65 

u.io) : u.ioo 

■ 'u.665 

12.651 


1.555 l.*^ 

1.8V ■ 1.7*0 
2.156 _ 2.015 
2.*6* 2.506 


U.665 
12.652 
i>*T 5 
I*. 525 
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, 2.815 

3.19a 

5.589 

*.012 

*.*60 

*.931 

5.12T 
! 5.9*6 

— 6.190 
7.057 

I 7.6*9 
8.26* 
6.90* 
.-9.57V 
10.265 
10.961 
U.725 
12.*95 
13.290 
1*.U2 


1.166 

1.739 

2.01I 

2.506 

2.621 

2.955 

5.510 

5.666 

1.065 

1.501 

1.91* 

5.105 

5.867 

6.590 


2.957 
5.515 
, 5.691 
1.090 
1.510 
*.952 
5.*i* 

_5-698 

6.105 

6.9J1 6.915 

7.177 ■ 7.*60 
8.017 : 8.027 
8.638 6.615 

9.251 9.25I 

9.885 j 9.85* 

-10.515,10.507 

U.221 • U.I6I 
'U.925 I U.877 
12.616 12.597 


0.067 
0.118 
; 0^155 
0.195 
0.258 
0.267 
0.511 

lS;IS 

0.679 
. 0.762 
: 0'.&9 
i 0.9S9 
■ 1.015 
1.100 
I 1.2W 
} 1.30 
: 1.12 
1.529 
, 1.652 

* J.79I 

1.951 


0.050 

O.OT5 

0.099 

• 0.129 
: 0.16I 

• O.805 

■ 0.217 

0.2S77 
0.?52. 
' O.lU 

I 0.160 

• 0.519 

0.585 


. 0.t772 
0.098 
0.128 
' 0.165 
' o.aoe 
0.215 

0.291 
; 0.519 
0.106 
! 0.151 

0.512 

0.571 


0.650 

0.721 To.708 


0.79* 

0.671 

0.955 

ifi- 

1.218 

1.559 

1.175 

1.599 

1.727 


0.779 

0.855 

0. 955 
1.022 

1. U7 
1.217 
1.325 
I.I3I 

1.552 

1.675 



0.51* 


o.p) 


0.399 


0.595 

. 

0.156 


0.152 


0.519 


0.510 

1 

0.581 


0.571 

.i„. 

0,6*7 

6.717 


.0.64^ 

. 0.70* 

1 

0.390 


0.775 

I 

0.666 


0.8*9 


0.9*6 


0.980 

- 

l.op 

1.122 


1.000 
’ i.096 


1.216 


).l86 


1.U7 


1.202 


l.*19 


' 1.P9 


1.526 


Kiai 


}.6p 


i- 1.507 


1.752 

l.lO* 

1.6y? 


i.inc 

).8o6 

i.floi 


1.«?96 

1.922 

i.920 

?,?)9 

?. ;?* 

?.0l? 

; i-.OfcO 

2-P7 

2-257 

2.166 

' 2.U»5 

2. 518 

2-395 

2.295 


C.67* 

2.530 

?.1?1 

r.kp 

?.B9i 

C.603 

2.558 

C.55* 

3.001 

2.831 


2.692 

3.i7? 

C.9«n 

c.sie 

f-.8‘3 

5.3*0 

3.751 

P.-/9* 


>.5?e 

>. 316 

3-153 

5. 127 

3.715 

J.107 

3-f«7 

3-279 

3. pi 

>.66? 

5-*U 

>.*55 

1.100 

3.0*5 

5-to5 

3-595 

*-301 

1.0?9 

3-769 


l.'/07 

*.XO 

i-'jyr 

s.«6 

*.7» 


l.U) 

i.<i90 

1.9P 

*.6ro 

1.286 

*.?72 

5.156 

1.620 

1.167 

1.151 

5.38* 

5.0*3 

1.551 

4.h55 

5.615 

5.:6? 

1.6k' 

*.0.73 

5.850 


5.0P 

5.011 

6.090 

5.7.r’ 


■...200 

6.1P 

5.'/X> 


5. 400 

6.507 

6.20b 

5.bp 



6.156 

5.0*6 

5.816 

7. 102 

6.7)1 

6.059 

6.027 

7.J68 

6.979 

6.275 

6.2*1 

7.6p 

7.?*9 

6.195 

6. *57 

7.912 

7.526 

6.72H 

6.60r> 

0. 700 

7.810 

6.9*8 


8.172 

6.103 

7.161 


8.750 

e.ioi 

. 7.11* 


.?•«*$ 

8.708 

7.65* 

7.60k 

9.3l« 

9>0?1 

7.098 

7.8*3 

9.6*9 

9. Pl 

8.1*5 

8.007 

9.959 

9.670 

8.396 


jlO.275 

10.006 

8.651 

, 6.58* 

10.396 

10.35* 

e.po 

6.6*0 

10.918 

10.706 

9.173 

9.098 

U.2*9 

U.069 

9.U1 

9.360 

U.562 

u.*» 

9.7U 

9.689 

U.925 

U.623 

9.907 

9.896 


12.210 

IO.S66 

io.no 

12.620 

12.6U 

10.5*9 

10.U6 

12.977 

15.021 

10.6P 


15. 5P 

13. Ul 

U.l^ 

U.01* i 

. )5t7P5 

iy07) 

U.*2* 

U.P*i 

ll.OTT 

11.516 

U.72* 

U.997i 

l*.*59 

Il.TTO 

12.029 

U,09*| 

i*.ep 

15.253 

12. 5^ 

12. IP 

15.827 

15-712 

12.650 

12. *99 

15.620 

16.200 

v.m 

u.aa0; 

16.022 

16.700 

13.291 

13.120 1 


JT-Rb 

38- Sr 

39- 7 

; lO-Zr 

ii-n 

i2.ito 

l>.Tc 

U.R, 

l5'Ri 

16- M 

17- i« 
id^d 

I9. tD 

50.80 

51- Sb 

52- Te 

55- 1 
5i-r* 
55<» 

56- 80 

57- U 

58- re 

59- Pr 

60- Rd 

61- ft 

62- aa 

63- EH 
61.cd 

65 - Tb 
ttrDy 
67-b 

66- Er 

69- -to 

70- Tb 

71- liJ 

I 72-8.' 
7VT0 
7I-W 
75-lte 

76-08 

TT-Ir 

78.pt 

79-08 

.60-% 

81-71 

62-1% 

Byn 

61. lb 
65-*t 
86 -Rb 
67-97 
8B-0i 
69.28 

98-0 

IK 

95-*i 

96-(k 


15.201 
16.106 
X7.0JT 
17.998 
. 18.967 

S).002 

21.051 

22. U8 
2J.221 

21.5*7 

25.517 

26.712 

27.980 

29.190 

30.186 

51.009 

35.16b 

31.579 

55.959 

37-110 

38.951 

40.1.19 

*1.998 

I3.571 

*5.207 

*6.6*6 

16.515 

50.229 

51.990 

. 53.789 
55.615 

57-W3 

59.382 

61.303 

63.30* 

65.313 

; 67.100 

69.506 

71.662 


15.161 
, 16.003 
1 17.011 

17.969 
1S.951 

19.961 

21.012 

22.072 

25.169 

‘26.6*1 • 
; 27.059 ■ 
29.106 ; 
3D.367 
! 51-690 
33.016 • 
:3*.398 
35.019 
37.255 ' 
’ 38.720 
*0.?3l 

*i.n2 

*5.3*9 

,**.955 

*6.561 

* 8 . 2*1 

*9.961 

'51.737 

'53-*9l 

155.292 

57.088 

50.969 

60.959 

62.9*6 

6*.9^ 

66.999 

69.090 ■ 

71.2WJ 


11.960 ' 15 
15.8>* • I* 
16.736 il* 
17.666 i 15 
16.621 I 16 
19.607 ! 17 
20.612 ' 16 
21.655 ; 19 

22.721 20 
23.816 ' 21 
2*.**2 22 
26.093 25 
27.27* 

a6-*85 .. 

29.723 26 
30.993 ri 
32.292 20 
33.6** ^ 
5*.96fc 
36.376 . 
37.799 55 
39.255 5* 
*0.7*6 
*2.269 
*5.8U ^ 

*5. *00 *0 
*7.027 11 
*6.718 12 

50.391 n 

52. IB? *5 
53.93* *7 
55.69? I9, 
57. *87 50. 
59-352 52 
61.262 5* 
65.309 55 
65.210 - 57. 
67.233 159- 
69.290 161 


1.39* 13. 3» 
>.16* l*.097 

.957 11.882 
i.n* 15.690 
>.61* 16.520 
.*70 17.373 
1-370 16.250 

L278 19.II9 
1.21* 20.072 
1.175 21.018 
‘.162 21.-. A 
1.172 22.W.' 

.207 21.000 

'.2T. 25.012 
••>57 26.109 
.171 27.200 

1.610 2«.ll5 
►.779 39.163 
'.97'’- p.623 
'.191 ».815 
i.Ho 33.053 
.717 5*. 276 

.023 55.518 
.559 36.815 
.726 38.180 
.12* 39.523 
.529 10.877 

.985 12.200 
.170 15.737 
.985 15.193 

.528 *6.686 
.099 18.375 
.730 *9.762 
.3b0 51.326 
.063 , 52.959 
.757 *5*.579 
521 156.270 

310 57.973 

IJl 1 59.707 


2.067 

2.221 

2.369 

2.517 

2.706 

2.881 

5.05* 

3.236 

3.119 

5.617 

3.810 

1.019 

4.237 
1.161 
1.697 
1.958 
5.190 
5-152 

5.720 

5.995 

6.263 

6.56) 

6.816 

7.JH 

7.1*8 

T.75I 

8.0^ 

8-393 

8.721 
9.003 
9.HI 
9.77« 

10. IW 
10.186 
10.867 
11.26* 
11.676 
12.090 
[12.522 


97-n 
9^ 
, 99-t 
lOO-M 


f 73.860 73.593 

I 76.097 75.605 

) 78. W9 T7.066 

I 80.713 80.165 
I 65. 106 82-2^ 
j^517"BT9W 
88.001- 87.2*5 
90.521 89.853 
95.U2 92.586 
95-T*0..9*.976 ; 

98.*li 97.626 ■ 
101.1*7 100.309 
] ioy9E?r 103.0*8 
' U6.759 105.038 • 

!j»-!asj08.6u . 

i U8.5ti;m.575 ! 
U5.59liU*.5*9 ' 
U8.619 U7.533 ' 
la.rao UD.5921 

! .«*>« 125,706 : 

; 128.0^ 126.875^ 
151.357 190.101: 
15*.685 135.^3 
! 1«.067 li.72* 

• 1*1.510 1*0.122] 

1 I 


7l.*0* 62.991 

75.5*9 6*.ea6 
75.756 66.820 
T7.960 66.79* 

80.98 TO.Cei. 

82.550~ 72.fio- 70I62 d| 
61.9B2; 7*.957j 72.79*' 
87.595 TT.097 7*.605 
89.609; 79.2961 76.8681 
92-519.' 81.525' 78.956j 
9*.877l 65.8001 Sl.oaol 
97.*85/ 86.U9; 85.6*5 
100.156 88.*85 ; 85.**6. 
102.816i 90.89*1 87.681. 

S .992 ; 9>.5J*i 89,9*2. 

.*08 95.W 92.?n' 
1U.S09' 98.*^ 9*.6ia' 
U*.lSl. 101.009! ^.025 
UT.1*6:105.655^ 99.*57, 
1^163 i08.351'101.952 
125.255 109.098 lOk.lld 
126 . 362 : m.896' 107 .025 
129.5 **!u*. 7*9 109.605 
152.76l;in.6*6 U2.2** 
15S.075il2D.596 U*.9e6 


12.965 

15.113 


15.073 


61. *77- 

65.270; 

65. Ul' .. 
66.980) 11.353 I 
68.6^! 11.611 I 

*-15.5*<I 
15.070 
16-993 
16.938 : 
17.I90 

18.058 ! 
28.658 ' 

19.255 i 

19.8*2 
20.160 i 
21.102 ‘ 
H.755 - 
22.117 ^ 
25.097 ! 

• ??-T95 .. 
rf.909 

25.250 
25.971 
26.729 

*7.505 


I 1.666 
! 2.006 
' 2.15* 

, a.505 
2.*67 
! 2.6CT 
2.795 
2.966 

5.1*5 

3.529 

3-520 

3-777 

5-939 

*.157 

1.381 

*.613 

1.656 

5.10* 

5.358 

5.625 

5.89* 

6.165 

6.k*3 

6.777 

7.OJ8 

7.312 

7.62* 

7.9*0 

8.258 

8.591 

e.920 

9.263 

9.628 

9.9T7 

10.^5 

10.73* 

11.130 

iU.555 

|U.955 

: 12.365 

12.619 
i 1V268 
15.753 
i*.ai2 
11.697 

15.607 

15.716 

I6.9I* 

16.78* 

17.357 

17.90* 

10.*8l 
19.078 
19.668 
80. 5U 
80.9*5 
21.596 
22.862 
22.9** 
25.6*0 
2*.552 
25.080 
25.621 
86.58* 


1.006 

1.9*1 

2.CT79 

2.220 

2.37* 

2.525 

2.6TT 

2.837 

3.002 

3.172 

5.352 

3.558 
>.729 
5.926 
1.132 
1.3*1 

1.559 

*.782 

5.OU 

5.2*7 

5.109 

5.72A 

5-960 

6.215 

6 M 6 

6.721 

6.983 

7.252 

7.519 

7.790 

8.07* 

6.361 

8.652 

8.9*3 

9.?ll 

9.556 ■ 
9.076 

10.198 
10.551 I 

■ 10.869 
; U.2U 

1 U.559 
U.919 
( 12.« 

• 12.657 
15-0** 

' 1V*2* 

I 15.817 

: 11.215 
> II.61B 
- 15.0tf 
IS. 1*2 
15.865 
16.296 
l6'.T3i 
27.165 

17.61* 

16.066 

18,525 

18.990 

19.161 

19.9^ 
20.122 
20. 912 


2,502 
2. *62 

2.6? 3 
2.792 
2.9f‘* 
3.1*4 
5.320 
3.5J9 
5.716 

>.920 

*.13) 

1.3*7 

1.570 

b.aoo 

5.036 

5.f«0 

5.551 

5.709 

6.05? 

6.322 

6.60c 

6.091 

7.100 

7. *78 
7.788 
8 . 10 * 

8. *tR 
8.7*8 

9.OS9 

9. *2l 
9.TT9 

10.112 

10.51* 

IO.B92 

11.285 


80.77* 16.837 
21.*01 17.25* 
22.0*2 17.677 
22.699 18.106 
25.370 18.5*0 
2*.056 18.980 
21.7^ 19.*26 
25.*75 19.r 


079 


17.TI0 

10.278 

18.629 

19.399 

19.970 

20.562 

21.166 

0.785 


1V615 15.*3S 
1V9*5 15.758 
11.979 11.089 
11.616 ll.lU 
11.961 1*.7*5 
15.909 15.079 
15.661 :i5.*20 
16.018 15.76* 
16.379 .16.U3 
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TYPICAL MAIN PHOTON ENERGY DISTRIBUTION MATRIX 





























PAIR PRODUCTION ENERGY DISTRIBUTION MATRIX 
FOR THE MAIN PHOTON FIRST COLLISION 
SHOWN IN FIGURE 9 
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PAIR PRODUCTION ENERGY DISTRIBUTION MATRIX 
FOR THE MAIN PHOTON SECOND COLLISION 
SHOWN IN FIGURE 9 
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TABLE 7 


ENERGY ABSORPTION TERMINATION DISTRIBUTION FOR 

FIGURE 9 HISTORY 



Photon 

Energy 


Weight 

(mc2) 


0.02 

8.0 


0.20 

4.0 

MAIN 

0.05 

4.0 

HISTORY 

0.15 

5.9 


0.20 

5.9 


0.10 

5.4 


0.07 

7.0 


0.03 

7.5 

FIRST 

0.09 

7.5 , 

PP 

0.03 

7.2 


0.04 

7.2 


0.02 

7.2 


0.05 

1.0 

SECOND 

0.06 

0.2 

PP 

0.07 

0.8 


0.01 

0.8 


0.01 

0.8 
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TYPICAL VALUES OF THE CHARACTERISTIC 
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TABLE 10 

COMPARISON OF MONTE CARLO AND ANALYTIC SINGLE SCATTER 
TOTAL NUMBER ALBEDOS 03 $ 


MATERIAL 

ENERGY 

(MeV) 

d 

(cm) 

L 

(cm) 

'“fes'aMl. 

(^ss) 

aj)'mc 

Xlff2 1 


.622 

■ 

17.0 

6.56 

6.42 

Al 


HM 

15 

5.47 

5.40 


1.25 


17.0 

4.71 

4.79 



■i 

4.25 

4.33 

4.20 


.662 


6.0 

4.29 

4.48 

Fe 


Hi 

1.19 

3.61 

3.22 


1.25 



5.82 

5.77 



■1 

1.5 

5.37 

5.61 


.662 

B 

5.0 

2.92 

2.89 

Sn 


Hi 

0.752 

157 

2.76 


1.25 


7.0 

2.48 

2.13 



B 

1.75 

2.38 

2.28 


.662 



1.12 

1.30 

Pb 


^H|H 


1.05 

0.97 


1.25 



1.02 

1.06 



H 

1.0 

1.00 

1.22 
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TABLE 11 


PHOTOFRACTIONS AND EFFICIENCIES FOR 
3" X 3" Nal (Tjt ) UNCLAD AND CLAD CRYSTALS 
EXPOSED TO A 10cm DISTANT POINT SOURCE * 

(1.18 g/cm2 ABSORBER) 

Photon Energy Photofraction Efficiency 


(MeV) 

Pu 

Pc 


^c 

4.45** 

.280 

(.278) 

— 

.727 

(.726) 

— 

2.754 

.268 

(.254) 

.251 

.423 

.434 

2.615 

.274 

(.262) 

.256 

.426 

(.428) 

.433 

1.276 

.402 

.365 

.506 

.519 

.835 

.500 

.449 

.576 

.584 

.6616 

.571 
(.562) 
(.569) ( 48 ) 

.507 

.615 

(.619)*^’ 

.622 

.511 

.662 

.588 

.661 

.667 

.32 

.845 

(.832) 

.749 

.766 

— 


* Unclad reference data indicated in parentheses. 

** 4" X 4" Crystal and Parallel Beam 
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Brief Description of the OSS Code 
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APPENDIX I 


/r 0\ 

Brief Description of the 058 Code' ' 

058 is a Monte Carlo code which directly simulates the experimental 
techniques used to obtain the pulse height distribution for a parallel beam 

/r 0\ 

of mono-energetic neutrons incident on organic scintillator systems^ ' 

The code is patterned after the far mare general 05R code^^^^ and, 
for the sake of brevity, at least passing familiarity with that code is 
assumed. Due to the nature of the specialization which 058 represents , 
a considerable amount of the memory storage required for 05R is avoided 
allowing the analysis of neutron events as they occur rather than after all 
have occurred. 

Although the code is capable of calculating the pulse height 
distribution for neutrons in the energy range 20 KeV to 70 MeV, experimental 
verification presently exists only for the range 200 KeV to 20 MeV. The 
code is written in the FORTRAN IV language for the IBM-360/91 computer. 

In this appendix some general features as given in reference (5) 
are repeated. An example input/output sample is given in Appendix V of 
report Volume II 
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INTRODUCTION 


In the calculation of the response function of an organic scin- 
tillator to neutrons of all energies it is necessary to obtain pulse height 
distributions for mono -energetic neutrons at energies where experimental 
determinations are often quite impractical because of problems with back- 
ground and absolute monitoring of neutron flux. For this reason it is 
desirable to have a method for accurately predicting the pulse height 
distribution per unit neutron flux over a wide range of energy and at 
reasonably small increments of energy. Because it allows direct simulation 
of typical calibration experiments with mono-energetic neutrons while 
utilizing all available empirical information (e.g., kinematics, angular 
distributions, cross-sections, light yields), a Monte Carlo "experiment" 
is frequently used for this purpose. 

The Monte Carlo "experiment" follows the course of each neutron 
through the scintillator and obtains the energy-deposits of the ions pro- 
duced by elastic scatterings and reactions. The light pulse produced by 
the neutron is obtained by summing up the contributions of the various 
ions with the use of appropriate light ys,. ion-energy tables. 

To a large extent the 053 code is patterned after the far more 
general neutron transport 05R code^^^^. The latter code covers the energy 

-3 

range from 77.13 MeV to 0.070 x 10 eV. The energy range is divided 
into 40 supergroups , each a factor of 2 apart; and each supergroup is 
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divided into 2^ subgroups, where n=0, 1, 9, and n is specified by 

the user. Since 05R allows the treatment of complex geometries consisting 
of as many as 16 media with any practical number of nuclides per medium, 
a tremendous amount of cross-section information may be required. The 
code stores in memory the scattering probabilities for one supergroup and 
processes neutrons through that supergroup in batches of from 200 to 
2000 neutrons. The scattering probabilities for the next lower super- 
group are then read into memory and the process is repeated. The product 
of an 05R run is a magnetic tape(s) which contains information describing 
each reaction which occurred. The results of the Monte Carlo experiment 
are found by analysis of this tape(s). 

Since the OSS code is concerned with only one medium, namely: 
the scintillating material, consisting of two elements and a total of eight 
reaction-types or scatterers, in the energy range 70 MeV to 200 KeV, the 
amount of cross-section information is small enough that the information 
for all supergroups may be stored in memory simultaneously. . The 
result of interest is that neutrons may be followed Individually, their 
histories analyzed as they occur, and, upon completion of the experiment, 
the results analyzed to obtain the pulse height distribution during one pass 
on the computer. 

The other fundamental differences between the two codes are that 
058 does allow the absorption of neutrons, but does not allow splitting or 
Russian roulette (biased weighting schemes). 
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A discussion of the use and experimental verification of the cal- 
culations of 058 is presented in a separate report on the calibration of the 
scintillator^^^^ . 

This IBM-360/91 version of the 058 code was obtained from the 

Radiation Shielding Information Center (R8IC), Oak Ridge National Lab- 

( 0 ) 

oratory' . Reference (5) consists of a 'user's manual' for the 058 code, 
the X8ECT code must be run to prepare the cross-section tape for 058. 

The X8ECT code is part of the 05R package obtainable from R8IC/ORNL. 

Both 058 and X8ECT were made operational for NA8A/GSFC by NU8 Corporation. 

A general treatment of the physics and kinematics of the organic 
scintillator as a neutron detector is given in reference (53). 
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OSS: General Features 


Energy Units and Cross-Sections 

Input and output energies for OSS are conyentional eV units. 

Internal calculations are made with energy represented in neutron speed- 
2 2 

squared units (cm /sec ). ' The conversion for neutron energy in eV to 

2 2 

speed squared in cm /sec is 

1 cm^/sec^ = 1.913220092 x 10^^ eV. 

OSS requires mean-free-flight paths and cumulative reaction 

( 51 ),^ 

probabilities in the manner provided by Code 6 of XSECT ttne cross- 

section handling code of the OSR package). Presently OSS uses the 12 

energy supergroups from 77.13 MeV (2 cm'^/sec'^) to 18.83 KeV 

55 2 2 

(2 cm /sec ) with 128 subgroups per supergroup. 

Scintillator System 

Two geometrical situations for cylindrical scintillators are 
presently available. The first as in Figure I-l Is brie of considerable 
interest . A parallel beam of mono-energetic nuetrons is 
incident on the portion of the lateral surface of the scintillator de- 
noted by the negative x-axis. 

In the second situation, as in Figure. 1-2 j . the beam_ of. neutrons is 
incident on a concentric circular portion of the flat surface of a 

scintillator system consisting of: - — 

(1) A Cj detector whose radius is that of the scintillator and 
whose thickness is given in equivalent, units, of the scintillating material 
This detector is often used in experiments to prevent the entrance of . 




SCINTILLATOR SYSTEM WITH NEUTRONS INCIDENT ON LATERAL SURFACE 


ORNL-BVro 67-12335 


z 



Figure I-l 
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SCINTILLATOR SYSTEM WITH NEUTRONS INCIDENT ON FLAT SURFACE 


ORSL-DTO 67-12336 



Figure I-2 
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charged particles into the scintillator and tends to slow down incident 

, (54r 

neutrons. 

(2) The scintillator proper. 

(3) A light pipe whose radius is that of the scintillating 

material. This unit has the effect of scattering neutrons back into the 

( 5 ^ 

scintillator. 

Generation of Neutron Histories 

A neutron history is generated in much the same manner as is used 

in 05R. 

(1) The neutron is given a NAME, a unique integer distinguishing 
it from other neutrons within the same case; initial energy 
and direction, and spatial coordinates of the point of 
incidence. 

(2) The path length to the next reaction point is chosen by 
multiplying the mean free patH for the current neutron 
energy, by an exponentially distributed random number. The 
coordinates of the next reaction point are then calculated. 
When this point is outside the system, the neutron has 
escaped and analysis of the neutron history inside the 
scintillator is performed. 

(3) The reaction which is to occur at the point chosen is selected 
by a random process from the possible reactions, given in 
TableJ"! according to the various partial cross-sections. 

(4) The appropriate reaction routine calculates the parameters 
associated with the neutron after reaction# unless absorption 
occurs# :and the energies of the reaction products. 
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TABLE I-l 


NEUTRON-HYDROGEN, NEUTRON-CARBON REACTIONS AND PARAMETERS 


Reaction 

Q-Value 

(MeV) 


Threshold 

(MeV) 

H(n,n)H 

0 


0 

12_, -12_ 
C(n,n) C 

0 


0 

12_. , ,12_ 
C(n,n y) C 

4.43 


4.8 

12 9 

C(n,a)^Be 

5.71 


6.18 

^^C(n,n'3a) 

12.73 

(56.93%) 

13.8 


11.1 

( 5.14%) 

12.01 


9.63 

(37.93%) 

10.42 


7.656 


8.3 

C(n,p) B 

12.613 


13.69 

12^, ai„ 

C(n,pn) B 

15.957 


17.29 

12 11 
C(n,2n)“c 

18.72 


20.13 


Note: Q-values (always negative) are used as positive in OSS. Since 

a pulse shape discriminator is used to eliminate y-pulses ex- 

12 12 

perimentally, the y from the C(n,n'y) C reaction is ignored. 
12 

For the C(n,n'3a) reaction the Q-values and branching ratios 

are those found to give best agreement with 14.4 MeV results. 

12 

The 7.656 MeV Q-value for C(n,n'3a) is used only in the rare 
cases where the cross-section is finite but the other Q-levels 
are energetically impossible. 



(5) The light yield of the reaction products is calculated unless 
the reaction occurred outside the scintillator proper but in- 
side the system. 

(6) If the neutron value has vanished (by absorption, escape 
from the system or loss of energy to less than 20 KeV) the 
neutron history is analyzed. Otherwise the neutron continues 
to be followed through the system. 

Reaction Analysis 

The light yield from each product of a reaction is calculated 
from the appropriate light vs. energy table. Presently there are tables 

for proton, alpha particle and nucleus light. The light for a ^Be, 

12 11 11 12 
B, B or C nucleus is assumed to be the same as for a C nucleus 

of the same energy. 

Since proton escape is considered two energies are calculated for 
the proton after the reaction - the energy with which the proton left 
the reaction and the amount of that energy not lost in the scintillator. 
The amount of light lost by the proton inside the scintillator is the 
difference of the light yields corresponding to the two energies. 

The light yield from each of a neutron's reactions is summed 
over all its reactions. • 

Neutron History Analysis 

For each neutron history the total amount of light produced 
during all its reactions is stored as a single pulse in a histogram. The 
light interval £ in which the pulse is stored is given by 
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i = 1 + A®‘^/0.015, 

where A is the total light yield and the constant 0.015 is called the bin 
width factor. The width of the interval ^ to in light units 
(cobalts) is given by 

B^ - ®£-i = (O.OISH)^ - (O.OlS(il-l))^. 

One cobalt is the base-line intercept of a straight line fit to the high 
end of a ^^Co gamma-ray pulse height distribution (linear plot). 

Pulse Height Distribution Calculation 

The histogram ("unsmeared" pulse height distribution) which re- 
sults from the Monte Carlo experiment is composed of light intervals 6A 
which vary proportionately to the square root of the light A. The final 
(smeared) pulse height distribution is composed of light intervals 6A in 
which 6A/A is constant. Each interval of the final pulse height distri- 
bution represents the integral over that interval of the initial histo- 
gram convoluted by a Gaussian function of specified width, normalized by 
the number of incident neutrons and the interval width. The result of 
this process is a smeared distribution which matches the smearing due to 
photomultiplier statistics, and is absolutely normalized to one incident 
neutron. 
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APPENDIX II 


The Kahn Monte Carlo Technique of Selection 

(l 8) 

of Scattered Photon Energy and Direction' 



22, Picking from the Klein-Niahina Scattering Formula 


It can be shewn that if a y r&y of energy a (in ec^) is 

traveling in a cloud of free electix>ns, then it has a probability 

g(x,a)dx<D. of having a collision in the diotance dl and emerging 

from this collision vlth an energy a* such that 

a . a a 1 

- < a» < ~ dx 

X- -X / 


r""V 


g(x,o)dx «" 


(co3^ 0-1+X+i)^ 

X 


1 < X < 1 + 2a 
X < 1, 1 + 2a < 


where cos 6 « l/a « x/a + 1 “ cosine of angle of ■ scattoring. 
a^ *» classical radius of electron 

- 2.81833 X 10“^^ cm 

“ ,2h9$}6 x cm^ 

n “ number of electrons per cc 

a’ » xa “ energy of y ray after collision. 


1-04- 


l*2a 


If a(a) 


g(a^x)dx then — 


is the so-called Congjton 


scattering cross section and is tabulated in R-170 Tables $ 

and 6, The probability that the y ray will have any kind of 
scattering collision in the distance dl is then a(a)dl, g(a,x) 
can now be written in the forms 

g(a,x) “ na(a)f(a,x) 

where f(cL,x) (cos^G - 1 + x + ~)A(a) 


-l+2a 


and 


K(a) 


(cos'^e - 1 ♦ X + 




dx 

vir 


f (a,x) is then the p.d.f . for the energy ratio x of a y ray whi.ch 
enters collison with energy a. The normalizing factor, K(a), is 
shown in Graph U« 

Graphs and Tables of the equation 


R 



f (a,x)dx 


can also be found in R-170. 

Engpirical fits to these functions could be made, but it is 
undesirable to have to fit a tv/o variable function unless absolutely 
necessary. 

Several exact methods for picking out of the Klein Nishina p.d.f, 
will now be discussed. This is partly psde.gogy and partly because 
the method used may depend on the particular section being picked 
from or the available machine. 


1 Actually cos G and l/a» as a function of R and a are tabulated, 
not X, 
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First, techniques in which x is drawn from a p.d.f. ra(x) 
and then accepted if R is < T(x) will be considered, TMs implies 
that E . f(a,x) » m(x)T(x), 


a. Let 


m(x) 

T(x) 


1 

(cos^0 


1 < X < 1 ♦ 2a. 


1 « X * i) 


It is clear that the expression to the right of the proportion 
2 

sign is < 2, since — < 0 and -^(x + i) < 2 in the i*egion 

“ x^ * 

1 < X < 1 ♦2a. 


So if 


T(x) 


^ (cos^0 - 1 + X ♦ ”) 


2x2 


the inequality 0 < T(x) < 1 is satisfied. The efficiency is 

„ „ m(x)T(x) 

- K(a) : 

“Ila 


^ • X/6^ 


m(x) 




1 < X < 1 4 - 2a 


T(x) « ^(cos^e - 1 + X ♦ i) 


y „ K(a) 

\ 5ihrr+2'^7 
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C. Let in(x) «■ ~ 1 < x < 1 ♦ 2a 

T(x) “ (coa^e - 1 + X ♦ i)/(l ♦ 2a ♦ 

I? - (l-«-2a)^ K(g) 

(2a) [l + (l*2a)^J 

That 0 < T(x) < 1 can be verified by noting that x ♦ J ^ 

ttonotonically increasing function in the region 1 < x and that 
2 

cos 6 ~ 1 has a maxiriuni value of 0. Therefore the function 
(cos^0 -l + x + i)<x+i<l + 2a + if 1 < x < 1 + 2a. 

d. It is also possible to break up f(x,a) into the sum of 
two p.d.f.'s; for example: 
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